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ABSTRACT 


Tnis paper examines the implications on acceptance 
sampling decisions when the Wald Sequential Probability 
Ratio (SPR) Sampling process is curtailed. Two procedures 
are proposed to determine the stopping rules. The first 
procedure uses the slope of the least-square fitted line 
compared with the slope of the boundary lines of a Wald SPR 
sampling Plan. The second procedure uses the relative 
position of the last observation between the rejection and 
acceptance lines to determine the stopping rules. Computer 
programs are used to simulate the sampling process, 


providing estimates of operating characteristic points. 
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เพ ร ON 


In a Wald Sequential Probability Ratio (SPR) Sampling 
Plan, sample size is a random variable and we can not 
determine the number of items to be inspected in advance. It 
could be a large or small number. This uncertainty could be 
prohibitive whenever the sampling budget is limited or time 
for a decision is constrained. In many cases it is 
preferable, for a variety of reasons, to have a finite upper 
bound on sample size. However, a general policy has not been 
found which determines how a final decision at the point of 
truncation should be made so as to conform to stated 
acceptable risks. 

The purpose of this paper is to examine the implications 
on acceptance sampling decisions when the Wald SPR sampling 
process is truncated by some predetermined sample number, 
and the final decision may therefore be based on statistics 
computed at the point of truncation. Two procedures are 
proposed to determine the decision rules for accepting a lot 
if sampling reaches the truncation line: (i) a least square 
fitted line method and (ii), a relative position of last 
observation method. 

In order to evaluate the implications of the proposed 
procedures on the risks associated with the sampling plan a 


computer simulation of the curtailed and uncurtailed Wald 





SPR sampling process is used, providing estimates of the 
probability of acceptance for various values of lot fraction 
defective. The Operating Characteristic (0.C.) curve of the 
curtailed and uncurtailed sampling are plotted in the same 
graph using & second Fortran computer program by the 
Versaplot-07 Plotting System available in Naval Postgraduate 
School Computer Center. 

The presentation starts with the nature of the problem 
which describes the method, the problem, and proposes two 
approaches to the problem. These are given in Chapter II. A 
description of the actual decision procedures and how the 
simulation was done is given in Chapter III. In the last 
chapter, the results of the simulation and the graphs of 


0.C. curves are discussed, and conclusions are drawn. 





TES NATURE OF THE PROBLEM 


In general, truncating a sequential sampling plan will 
increase the probability of type I and type II errors. The 
exact functional relationship between the size of error and 
the sample size of truncated sequential sampling is not yet 
known. However, its upper bound may be derived [B. K. Ghosh, 
ee. +, pp. 223}. 

The purpose of this chapter is to: describe the general 
concepts of Wald SPR Sampling Plan and its Average Sample 
Number; discuss considerations in the curtailment of 
sequential sampling; and describe two proposed procedures to 
determine decision rules for truncated sequential sampling 


plans. 


A. WALD SEQUENTIAL PROBABILITY RATIO SAMPLING PLAN 

Abraham Wald [ Ref. 9] simplifies tne process of sequential 
sampling by a scoring method with acceptance and rejection 
boundaries which will meet some preassigned requirements. 
If the score at any time becomes larger than the first 
boundary (i.e., rejection line) the lot is rejected. If it 
falls below the second boundary (i.e., acceptance line) the 
lot is accepted. There are four specification requirements 
which completely determine Wald SPR Sampling Plan for frac- 


tion defective. Those specification requirements are: 
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lc "E the acceptable quality level for the lot, 
expressed as a fraction defective, 

າ Pos the lot tolerance fraction defective, expressed 
as a fraction defective where "2 ን ງ່! 

3. «£ , probability of rejecting lots of quality ງງ " , 

4. ff , probability of accepting lots of quality po'. 

Graphically, a Wald SPR sampling procedure can be 
described as follows. Consider a chart which consists of a 
vertical axis representing the number of defectives, a 
horizontal axis representing the number of items inspected 
and a pair of parallel straight lines with positive slope 
which are uniquely determined by the specification require- 
ments. During the sequence of inspection the total number 
of defectives is plotted against the total number of items 
inspected on the chart. As long as the plotted points fall 
between two lines, the inspection continues. An inspection 
terminates when a plotted point falls on or outside either 
of the LP 

Defining upper line by R and lower line by A, where R 
and A are functions of sample number, the eguations of the 


lines may be written as 


E = no + sn 


and 


ሙ 
ዘ 


1 





where R will give a rejection number and A will give an 
acceptance number at sample number n. The constants S, hj 
and ho are the slope and the intercepts and their equations 


may be written as follows (Ref. ODD. 2.14] 


(1 -p ) 

۱۳ ... 
Em uS , (3) 
lbs) 
wr مس نت‎ 
๒ 0 ฝ 0 ไ ไ ฝ ว 
) 1 - A) 

108 
ภา = “ , (4) 
pi (1 - p) 
log 
PP 
and 
: (1- B) 
08 
A 
he © (5) 
p (1-p) 
<<< ኢ: 


ps ( 1 - pi ) 
In these equations, it is necessary that Dj must be less 


than p3 and <+ ñ is less than unity, so that quantities 


obtained before applying algorithms are always positive. 


B. AVERAGE SAMPLE NUMBER (ASN) 
Since sample size is a random variable, it is not 
possible to determine exactly how many items from a lot have 


to be inspected, but it is possible to compute the average 


a 


Mè 





depends on quantities h4 n, , and 8. The equations are 


as follows [ Ref. - ESE]: 


= ۳ (ny + 12) - ho 
p a ነ 
s- p 


n 


where P is the probability of accepting a lot of duality p', 


and ከገ ; h, and Ss are Computed from the specification 
requirements. In particular when p' - Dj we have 
1-&)n - & h, 
n ı = 9 
Pi s- p- ' 
when p' = Po , we have 
1-P)h,- Pa 
n << , 
፲ . _ 
2 Po S 
and when p' = S , we have 
(ny : ከ3) 
แพ cog QU ONE 
Note that pj' < 5 ረ ፔ2' and in general "py ረ ຖ : We 


normally observe an increasing average amount of inspection 
as p' goes from zero to P  » and a decreasing amount of 
inspection as p' goes from p5' to unity. Hence the greatest 
ASN is required for a lot with quality between ນງ and 2 ۰ 
In addition, the greater the risk sizes «ď and fò are, the 
smaller also the ASN. These properties are useful when we 


discuss the curtailment of sequential sampling. 
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Figure 1 shows a typical ASN curve. The vertical axis 
represents np , the average sample number, and the 


rizontal axis represents p', fraction defective of lot. 


Ro 
p 
(average 
sample number) 
i 





t : 
: - --p(fraction 
9 P. S Pa 1 defective) 


FIGURE 1. A TYPICAL ASN CURVE 
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C. CURTAILMENT OF SEQUENTIAL SAMPLING 

As mentioned before, the number of samples required in 
Sequential sampling to achieve a conclusive decision is not 
a constant but a random variable. In practice it may be 
desirable to have an upper bound on sample size. Setting the 
sample number constant could increase risk size, since the 
Sequential process may or may not terminate before the 
truncation point. Two steps could be considered: (i) Try 
to reduce ASN, (ii) Modify the sequential sampling plan. The 
first step may require a compromise among the quantities A , 
[9 > Pi , and Pa . À property of the ASN states that the 
greater KA and ይ or the larger the difference between ມງ 
and Do. , the smaller the ASN will be. We could then make 
adjustments either in the quality limit (i.e., pi and 22” 
or in the size of risks (i.e. ຫຼ and f ) or both. But this 
compromise might not be applicable if the specification 
requirements are strictly kept. The second step is suggested 
by J. J. Bussgang and M. B. Marcus in Reference l. Instead 
of using a pair of straight lines as the boundaries, they 
propose "gently sloping" lines as the boundaries so that 
they would monotonically converse as the sample number 
increases. 

In both steps the given point of truncation must be 
determined sufficiently beyond ASN so that most of the 
sampling terminates before the truncation point is reached. 
This is clear because otherwise probabilities of the first 


and second kind of errors will increase. The two above 


JE 





procedures have limitations on their applicability. Now, let 
us develop two simple procedures. 

The first procedure uses the slope of a least-square 
fitted line as an estimator for the direction of the plotted 
sequential sampling, and this slope is compared with the 
slope of boundary lines computed from specification require- 
ments. The second procedure uses the fact that when sample 
Size is sufficiently large, the total defectives will either 
close to the rejection line or close to the acceptance line, 
with the probability of eventually crossing either line 
equal to unity. Further discussion of the proposed procedures 


is given in the next section. 


| PROPOSED PROCEDURES 
Consider a Wald SPR sampling plan with specified Ko, 


(3 > ນ່ , and Do. . Let n' denote the maximum sample 
number to be allowed, which is determined before the 
sampling begins. 

1. Least Square Fitted Line Method 
suppose up to n' there is no decision made either to 


accept or to reject the lot. By then we have observed: 


(1,X,), (2,X,), (3X), '. ። ea (ASK) a » » › (ቪከ',ጄ ") 
Where n is the sample number and x is number of defectives 
found in n observations. A least-square line fitted from 
the origin through the observed samples will have a slope 


given by: 
16 





8 > ร 
n' 
n 
ก L 
Di 
where A, is number of defectives at stage n, and = n 
Ml 
can be simplified by 
n' : 
Poe | 
2 
PS 
wielding 
DE 1۳ tle m X 6] 
E 2 »5 ຈ” | 
n=1 


Let us compare the slope of the least-square fitted line b 
with the slope of the boundary lines s. The decision rules 

| 
are given as follows. If b is greater than s we reject the 
lot and terminate sampling. If b is equal to or less than s 
we accept the lot and terminate sampling. 

2. Relative Position of Last Observation Method 

Again, suppose up to n' observations no decision can 

be made. This means that Xa is always between the boundary 


We forn= 1, 2, 3, +... n'. For a lot with quality 


better than Dj» the number oí defectives tends to close to 


1% 





the acceptance line if n is getting larger. On the otherhand 
Mar a lot with quality worse than Do's the number of defect- 
ives tends to close to the rejection line if n is getting 
larger. Let us take a constant distance above acceptance 
line, denoted by d. We can then define a new acceptance 
number A', where 


; [?] 


The decision rules are given as follows. If At is greater 
than A', reject the lot and terminate sampling. If ແ is 
equal or less than A', accept the lot and terminate sampling. 
Implementation of either of these procedures will have an 
impact on acceptance probability, and the curtailed plan 
should have a different O.C. curve from the original uncur- 
tailed plan. The magnitude of the change of the O.C. curve 
may be evaluated by simulation. The simulation procedures 


are described in the next chapter. 


o 





III. EXPERIMENTAL PROCEDURES 


There are three distinct steps in acceptance sequential 
Sampling by attributes. First, determination of objectives 
Or specifications, second classification or good or bad 
items, and third a valid procedure of inspection. 

The experimental procedures discussed in this chapter 
are presented in accordance with those three steps, and then 
used to evaluate the implications of the two proposed 
stopping rules on the plan's operating characteristic curve 


| 111721 ท ๕ 8 computer simulation. 


۳ PROCEDURE I 

HT nding a Wald SPR sampling plan for fraction 
defective, the specification requirements d, ምዎ ዔዝ p,' , and 
Do are used to compute s, ከገ › and h. using Equations (3), 
(4) and (5), and these give the equations of the acceptance 
and rejection lines as functions of n. Now, we consider lots 
En auality p'. We draw items from the lot, one at a time, 
and classify each as good or defective, defining as the 
number of defectives found up through the first n items. If 
a is equal to or greater than the rejection number we 
terminate the sampling and reject the lot. If X is equal to 
or less than the acceptance number we terminate sampling and 
accept the lot. Otherwise we repeat sampling until n = n', 
where n' is the curtailment point. At stage n' we compute 


the slope b by Equation (6) and compare it with s. If b is 
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Emeater than s we reject the lot; if D is equal to or less 
than s we accept the lot. In both cases the sampling process 
terminates. 

In simulation of ihe use of this stopping rule, the 
overall sequential sampling is replicated through k lots, 
where k is a large number. The probability of acceptance of 
a lot of quality p' is estimated by the number of accepted 
lots divided by k. If we repeat the whole process for 
different p' then we will obtain additional points of the 
0.C. curve for this curtailed sequential sampling plan. 

The truncation points are computed before the sampling 
begins. In this paper n' is computed as percentages of nos 


since it represents the largest ASN. 


E PROCEDURE ILI 

The Procedure II is similar to Procedure I except that 
the truncation point, Kou is compared to an acceptance 
number A', where A' is a function of d and computed using 
1 on (7). By trial and error it turns out that for 
large n', the value of A' is approximately equal to sn' since 
d is approximately equal to hy. However, for small n' it 
gives poor O.C. curve. The stopping rules are then: if = 
is equal or less than sn', accept the lot and terminate 
sampling. Otherwise reject the lot and terminate sampling. 

The simulation of the two procedures at different 
truncation points was done simultaneously with the 
uncurtailed sampling. The details of the simulation are 


given in the next section. 
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C. COMPUTER SIMULATION 

Monte Carlo simulation was used to simulate the Wald SPR 
sampling process. The computer programs were written in 
Fortran IV and utilized the IBM-360 computer at the Naval 
Postgraduate School Computer Center in the period of October 
1979 to March 1980. 

Input variables consist of the four specification 
requirements (denoted by A, B, Pl and P2), the number of 
replications, and the number of points on O.C. curve. A 
uniform random generator (GGUBS) with double precision was 
used to classify as good or as bad an item from a lot. To 
save computer time, the simulation of both procedures each 
with 5 different truncation points and the simulation of the 
uncurtailed Wald SPR sampling were done simultaneously in 
one run for each pair of values for Pl and P2. Eighteen 
operating characteristic points were computed for each pair 
of values for Pl and P2, where Pl was given from one 
percent to ten percent and P2 was from five percent to 
thirty percent. The parameter values used to investigate 
the performance of each procedure are shown in Table l. 

A second computer program was written in Fortran IV to 
plot the 0.C. curve of uncurtailed and curtailed sampling 
in one graph, where the data points were obtained from the 
first computer outputs. This will provide a visual repre- 
sentation of the difference between the two 0.C. curves. 

Ihe plots were done by Versatec-07 Plotting System available 


in the Naval Postgraduate School Computer Center. 
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Michael W. Gavlak ( Ref. ພາກກ, 24-26 | stated that to 
simulate estimates of O.C. points for repeated Bernoulli 
trials with p' ranging from one percent to thirty percent, 
it is sufficient to take 5000 replications of each estimate 
within reasonable accuracy, namely two or three decimal 


places. 


TABLE 1. PARAMETER VALUES USED IN SIMULATION 


EMOD. of type I error pemos 

EMOD. OL type II error በን 

Acceptable quality levels cO. 0:295, 0 10 

Lot quality tolerance values . 00 10.1 500 30 
Number of replications : 0 

Number of O.C. curve points Ze 

memeentages of n.for curtailment 12  . 1 0 
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IV. RESULTS AND CONCLUSIONS 


In general when Wald SPR sampling process is truncated 
with the same stopping rule, then its O.C. curve varies as 
tae point of truncation varies. The larger the point of 
1 Lion, the closer its 0.C. curve to the O.C. curve of 
uncurtailed sampling process. Using no» the average sample 
number when lot fraction defective is s, as standard for 
comparison, the graphs show that for n' greater than 150 
percent of Ng» their O.C. curve gives good approximation 
to the O.C. curves of uncurtailed sampling, since most of 
the samplings terminate before n'. 

Comparing the results of Procedure I and Procedure II, 
the numerical output shows that for large n', Procedure II 
mies better approximation to the uncurtailed O.C. curve. 
1 ከ ከ ከይን” Procedure II is a more simple method, hence it is 
ore practical. The determination of constant d, however, 
needs further investigation, particularly for small n'. 

For further investigation, notice that the Procedure I 
which requires the least-square fitted line through the 
orlgin raises question whether an ordinary least-square 
fitted line will give better approximation even though it 
may be less practical. Another area for further study may 
include the possibility of using the variance of ASN to 


determine the proper location of. the truncation point. 
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In all, sixty cases were examined and twenty two of their 
O.C. curves were graphed. The tables of the first computer 
output and the results of the second computer program are 


presented below. 
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CONTINUE 
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CONTINUE 
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ห บ เห ค ว ร (16F29/316ZyóáFS (3 6 F5S,O 12) 


STOF 
END 
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